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Abstract The major sources of the Soft X-ray Background (SXRB), besides distinct struc- 
tures as supernovae and superbubbles (e.g. Loop I), are: (i) an absorbed extragalactic emis- 
sion following a power law, (ii) an absorbed thermal component (~2xl0 6 K) from the galac- 
tic disk and halo, (iii) an unabsorbed thermal component, supposedly at 10 6 K, attributed 
to the Local Bubble and (iv) the very recently identified unabsorbed Solar Wind Charge- 
eXchange (SWCX) emission from the heliosphere and the geocorona. 

We study the SWCX heliospheric component and its contribution to observed data. In 
a first part, we apply a SWCX heliospheric simulation to model the oxygen lines (3/4 keV) 
local intensities during shadowing observations of the MBM 12 molecular cloud and a dense 
filament in the south galactic hemisphere with Chandra, XMM-Newton, and Suzaku tele- 
scopes. In a second part, we present a preliminary comparison of SWCX model results with 
ROSAT and Wisconsin surveys data in the 1/4 keV band. 

We conclude that, in the 3/4 keV band, the total local intensity is entirely heliospheric, 
while in the 1/4 keV band, the heliospheric component seems to contribute significantly to 
the local SXRB intensity and has potentially a strong influence on the interpretation of the 
ROSAT and Wisconsin surveys data in terms of Local Bubble hot gas temperature. 
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1 Introduction 

The surprising discovery of X-ray emission from comets dLisse etal.Ul996h led to estab- 
lishing Charge-eXchange (CX) collisions between highly charged Solar Wind ions and solar 
system neutrals as a very efficient non-thermal mechanism of soft X-ray emission dCravensl 
ll997h . lCo3 Jl998h suggested that X-ray emission induced in Solar Wind Charge-eXchange 
(SWCX) collisions with interstellar (IS) neutrals flowing across the Heliosphere would be 
an additional component of the diffuse soft X-ray background (SXRB). 

Until that discovery, the SXRB was generally assumed to consist of an absorbed extra- 
galactic emission following a power law, an absorbed thermal component (~2xl0 6 K) asso- 
ciated with the galactic disk and halo, and an unabsorbed thermal component, supposedly at 
10 6 K , attributed to the Local Bubble dMcCammon and SandersL l 1990; Kuntz and Snowdenl 
l2000l . and references within). The Local Bubble (or Local Cavity) is defined as the region 
within ~100 pc of the solar system partic ularly devoid of dense ( and cold) gas as suggested 
by absorption lines in nearby star spectra (Lallemen t et alll2003l) . 

The 1/4 keV background generally presents a negative correlation with the column den- 
sity of Galactic neutral hydrogen (dominated by a Galactic plane to-pole variation) that 
strongly suggests an "absorption" model, in which the sp atial structure of the SXRB is pro- 
duced by absorption of a distant X-ray component (e.g., iBowver et allll968l : lMarshall and ClarkL 
Il984t) . However, the strong absorption predicted in such models was not confirmed by ob- 
servation, neither was the large energy dependence of the absorption, so different solutions 
were needed. 

At the time, the model most consistent with th e data was the so-called displacement 
model dSanders et all 1 19771 ; ISnowden et all [l990a). In this model, the hot gas is in the 
foreground (in the Local Hot Bubble) with respect to the absorbing cool-gas regions, and 
higher X-ray intensities are produced in directions where the cavity has a greater extent 
(towards the LB chimneys at high galactic latitudes) and thus, a greater emission measure 
of the plasma in it. 

Since the discovery of SWCX emission, it became clear very quickly that it should 
contaminate at an unknown level all soft X-ray observations and it should depend strongly 
on SW flux and abundance variations as confirmed by the Long Term Enhancements de- 
tected dur ing the ROSAT all-sky survey (Snow d en et all 1 1995b and recent X-ray obser- 
vations dSnowden etail |2004 ISmith et allT2005l ; iHenlev and Shelter! l2007al) . The debate 
concerned the amount of contamination, and how the SWCX emission level would interfere 
with the Local Bubble thermal emission, since the two are most probably the major com- 
ponents of the unabsorbed fraction of the SXRB. Post-CX geocoronal emission can also 
contribute, bu t only during very in tense SW enhancements and for short and easily identi- 
fied intervals dCravens et alll200ll) . 

The two components, although due to completely different mechanisms (collisional ex- 
citation for thermal plasma emission and electron capture for CX emiss ion), produce sp ectra 
in the same energy range (E < 1.5 keV) that are frequently confused. ICraven s (2000) esti- 
mated that SWCX emission could be as much as the LB emission. Line ratios, though, are 
very sensitive to gas temperature (and abundance) for thermal emission. For example, a 10 6 
keV plasma might be responsible for all IS diffuse emission at 0.25 keV, whil e at 2-4 x 
10 6 K it might significantly contribute to the 3/4 keV diffuse emission as well ( Raymond, 
Il988l ; IMcCammon and SandersL Il990l and references therein). CX line ratios mainly de- 
pend on heavy SW ion rel a tive a b undances and collisio n energy (Beiersdo rfer et all 1200 ll ; 
iKharchenko and Dalgarnol l200ll ; IWargelin et all [2008). Detailed spectral comparison of 
SWCX and thermal emission needs a spectral resolution that should allow to resolve char- 
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acteristic features as the O vn triplet at 0.57 keV and in general individual lines from the 
continuum. Such a comparison is beyond the scope of this paper. 

In this paper we mainly focus on the intensity contribution of the SWCX emission to 
the observed (foreground) data in the 0.1-1.0 keV range. In section|2]we briefly present the 
basic simulation model that calculates the SWCX emission. 

In section[3]we analyse the simulation results for the 3/4 keV range and in particular the 
oxygen lines (O vn at 0.57 keV and O vm at 0.65 keV). We present typical stationary maps 
of oxygen lines in section [3~T1 and in section [3^2l we present the dynamic model we developed 
to account for short-scale temporal variations of SWCX emission , with an application on 
the MBM 12 Suzaku observation of February, 2006. Finaly, we conclude section [3] with a 
summary of shadowing observations on which we applied our SWCX dynamic model and 
discuss the results in the 3/4 keV range (sect. l3.3l >. 

I section|4]we present a preliminary analysis of SWCX simulation results in the 1/4 keV 
range that we compare to observation data from the Wisconsin sounding rocket and ROSAT 
satellite surveys. 

Finally, in section[3]we close our paper with a general discussion on our results. 



2 SWCX model description 

In iKoutroumpa et alj d2006h we have presented the basic stationary model calculating the 
SWCX emission in the inner heliosphere. We calculate self-consistently the neutral H and 
He density distributions in the inner heliosphere (up to ~100 AU), in response to solar grav- 
ity, radiation pressure and anisotropic ionization processes for the two neutral species. Ion- 
ization is mainly due to charge-exchange with SW protons for H atoms and to solar EUV 
photons for He atoms. We also consider the impact of CX on the solar wind ions distribu- 
tions. This interaction is described in the following reaction: 

X Q+ + [H,He] -^X* (Q - l)+ + [H + ,He + ] (1) 

The collision rate per volume unit R x e+ (cm -3 s~')ofX^ + ions with the neutral heliospheric 
atoms is given by the equation: 

Rxa*(r) = N xQ+ (r)v r (o- (H xQ+) n H (r) + cr (HeX Q +) n He (r)) = R {X Q+ ,H)( r ) + R (xQ+ ,He)(. r ) ( 2 ) 

where cr^ HX Q+) and cr^ HeX Q+- ) are the hydrogen and helium CX cross-sections, n#(r) and 
nn c (r) are the hydrogen and helium density distributions respectively, v r = Vsw -v„ ~ Vsw 
the relative velocity between S W ions and IS neutrals in the inner heliosphere, and N x q+ (r) 
is the self-consistent solution to the differential equation: 

dN X Q+ 

dx = -N x a + {o- (HxQ+) n H (x)+ a- {HeXQ+) n He (x)) (3) 
+ A^(e+n+ (o- (HX (Q+i) +) n H (x) + o- {He X (Q+i) +) n He {x)) 

expressing the evolution of the density distribution of ion X^ + along SW streamlines due to 
production (from CX reactions of ion X^ +I >+) and loss terms. 

Cross-section uncertainties are mainly due to instrumental systematic errors and most 
important to collision energy dependance of cross-sections. Detailed uncertainties for indi- 
vidual ions are not giv en in literature, but average uncertainties of ~30% at most are reported 
dWargelinetaill2008h . 
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Fig. 1 SWCX model spectra in Line Units (photons cm 2 s 1 sr 1 ) for the 3/4 keV band. The emitting ion is 
noted above each line. 



Then, we establish emissivity grids in units of (photons cm 3 s ): 

Bi(r) =R {X Q+,H)( r ) Y (E i ,H)+R { xQ+,He)( r ') Y (E i ,He) (4) 

where Y(E h M) is the photon emission probability of spectral line E, following CX with the 
corresponding neutral species M (H or He individually). For any line of sight (LOS) and 
observation date, this spectral line is given by: 



I El (LU). 



1 

4^ 



I 



AOQAU 



Bi(s)ds 



(5) 



which defines the average level emission of the spectral line for the particular date and LOS, 
as well as the solar cycle phase (minimum or maximum) corresponding at this date. 



3 SWCX in the 3/4 keV Band 

An example of calculated SWCX emission spectra in the 3/4 keV (0.5-1.0 keV) range is 
presented in figureQ] The major emission lines contributing in this energy range are: 

(i) the He-like O vn multiplet, consisting of the following transitions: 

- the 2 3 Si triplet at 560.9 eV, called forbidden (06f) 

- the 2 3 P[ triplet at 568.5 eV, called intercombination line (06i) 

- the singlet state 2'Pi at 574 eV, usually refered to as resonance line (06r), and 

(ii) the H-like O vm Lyman-o- line at 653.1 eV. 

In this list we can also add the He-like Ne ix multiplet (905.1 eV, 914.7 eV and 922.1 eV), 
but actual observing instruments do not always allow the detection of this line. 
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3.1 Stationary Model Maps 



In figure |2] we present an example of calculated full-sky maps of heliospheric emission for 
line vm (0.65 keV) in ecliptic coordinates, on December, 5 (when the Earth is located 
downwind) of a typical solar maximum (upper panel) and solar minimum (lower panel) 



the 



<t»B 1.0 

EARTH DW «■ ■ ' = 




MAX 



year. Colorscales represent O vm line intensity given in units of 10 erg cm sr 
red colour corresponding to minimum and the blue to maximum values. We have removed 
from every map a data portion of 20° x 20° around the solar disk where no instrument can 
observe. 

What we need to point out here are 
the striking differences between maximum 

and minimum solar conditions. At low lati- 

tudes, near the solar equator where slow SW 
dominates, X-ray emission is more intense 
for solar minimum than for solar maximum 
because the neutral atom content is higher 
during solar minimum. During solar mini- 
mum, photo-ionization is less efficient and H 
and He are less readily destroyed by photo- 
ionization. Since gravitational pressure ex- 
ceeds radiation pressure, neutral H atoms 
have incoming trajectories that fill the ion- 
ization cavity. 

On the other hand, at high ecliptic lat- 
itudes, LOS's are mostly affected by fast 
wind during solar minimum, and the X-ray 
emission is dominated by differences of ion 
relative abundances with respect to maxi- 
mum cycle phase. Indeed, 8+ is completely 
absent from the fast SW and the O vm line 
is expected to be much fainter at solar min- 
imum than solar maximum, as it is seen in 
the maps in fig.|2] 

O has the same trend in relative abun- 
dance variations as 8+ between solar mini- 
mum and solar maximum. Nevertheless, in- 
stead of being completely absent in the fast 
wind, it is only strongly depleted ([0 7+ / O] 
= 0.03) with respect to slow wind ([0 7+ / O] 
= 0.2). Therefore, we expect that O vn line 
(0.57 keV) variations from solar maximum 
to solar minimum will also have the same 
trend as O vm shown in the maps. 

Two excellent examples of such differences between solar minimum and maximum 
conditions for the same LOS, are the MBM 12 and S outh Galactic Filament (SGF) shad - 
owing observations , analys ed by ISmith et al] ( I20051 . 120071) and by Henley etal. (2007); 
iHenlev and Sheltonl fc007bll3) respectively. These fields were thoroughly analysed for their 
SWCX contamination in Koutroum pa et aLl d2007) and the main outline and conclusions are 
presented in section [331 
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Fig. 2 Calculated full sky maps in ecliptic coor- 
dinates of the 0.65 keV O VIII line for solar max- 
imum (upper panel) and solar minimum (lower 
panel) conditions. The observer is situated down- 
wind (ecliptic longitude 75°). Colorscales are in 
units of 10~ 9 erg cirT 2 s sr . 
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The MBM 12 LOS points toward galactic coordinates (159.2°, -34.47°), which translates 
to (47°, 3°) in helioecliptic coordinates, very close to the ecliptic plane and inside the limits 
of the equatorial slow wind zone. Therefore, during solar minimum the MBM 12 exposures 
(Suzaku) will yield higher SWCX contamination than solar maximum exposures (Chandra) 
in the 3/4 keV range. 

In the case of SGF, at very high southern ecliptic latitude (353°, -73°) the effect is the 
opposite compared to MBM 12, because it is dominated by differences in the solar ions 
relative abundances. Therefore, the SWCX oxygen emission in the 3/4 keV range for the 
SGF field is expected to be much fainter at solar minimum (Suzaku observations) than solar 
maximum (XMM-Newton observations). 



3.2 Dynamic Simulations 

In lKoutroumpa et all j2007l) we presented a dynamic variant of the basic stationary SWCX 
model, used to simulate the heliospheric X-ray modulation due to SW proton flux and abun- 
dance variations. For each field simulated we use real-time flux and abundance measure- 
ments from in situ solar wind instruments (Wind, ACE-SWEPAM) and we model SW en- 
hancements as simplified Corotating Interaction Regions (CIR) of a step function form. The 
total width of the CIR is defined such that the total duration of the step function is the same 
as the measured enhancement in SW instruments. The step function amplitude is defined so 
that its total area is equal to the integral of the measured flux during the same period of time. 

At each instant we define the form of the CIR (as a Parker spiral) and its propagation in 
the interplanetary space, taking into account solid solar rotation (27-day period), the radial 
propagation speed, the "ignition" time on the solar disk towards each radial direction, and 
the total width of the CIR. Only in the cases of CMEs, we neglect the solar rotation, since 
radial propagation is dominating the CME structure. Depending on the CIR's width and 
propagation instant, we calculate the fragment of the LOS affected by the CIR and the local 
emissivity ej(r) modified by the SW proton flux and heavy ion abundance variations. We 
can then reproduce the temporal variation of the X-ray intensity levels during the periods of 
observation in simulated lightcurves for each of our targets. Abundance variations, can be 
correlated or anticorrelated with proton fluxes, so they can either emphasize or compensate 
for the influence of the SW proton flux enhancements. 



3.2.1 The Case of the MBM 12 Suzaku Observation of February 2006 

The observation of MBM 12 with Suzaku on February 3-8, 2006 was performed in two 
consecutive exposures, ON-CLOUD (3-6/02/2006, for a tot al of 23 1 ks ) and i mmediately 



consecutive exposures, ON-CLOUD (3-6/02/2006, for a tot al of 231 ks) 
after that, OFF-CLOUD (6-8/02/2006, for a total of 168 ks) dSmith et all 



2007). A detailed 



description of the simulation we performed on this field is given in lKoutroumpa et al. I ( 120071) . 

As we already explained, the MBM 12 LOS is located at (47°, 3°) of helioecliptic co- 
ordinates. In the February Suzaku observation this geometry is pointing directly toward the 
He cone, which is only ~2 AU away from the observer, where the emissivity on the LOS i s 
maximum. The observation geometry was presented in figure 1 o flKoutroumpa et al] ( l2007h . 
The cone, with a denser distribution of helium atoms, is acting as an amplifier to any SW 
flux variation. 

The MBM 12 Suzaku observation was influenced by a short SW perturbation recorded 
at the end of the ON-CLOUD pointing. We model this enhancement as a step function in 
particle flux of 5.2 x 10 8 cirT 2 s _1 propagating at a low speed (Vsw = 350 km s~') for 0.75 
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Fig. 3 Top: O VII and O VIII line 
simulated lightcurves (in LU) for 
the MBM 12 ON and OFF obser- 
vations with Suzaku during the 
period 03-08/02/2006. Plain line 
stands for the O VII emission and 
dashed line stands for the O VIII 
emission. Dots represent the mea- 
sured ON and OFF O VII line flux, 
and triangles represent the mea- 
sured ON and OFF O VIII line flux 
temith et alll2007h . 
Bottom: Solar wind proton flux 
(dotted line) in units of 10 s cm~ 2 
s for the same period. The step 
function simulating the SW en- 
hancement is presented by the 
plain black line. In both panels, the 
vertical plain lines represent the 
start and end of the observation pe- 
riod, while the dashed vertical line 
is the separation between the ON 
and OFF exposures. 



05/02/2006 06/02/2006 07/02/2006 08/02/2006 

Date 



d. The SW proton flux remained at a high, but stable, level equal to 4.42 x 10 8 cirT 2 s _1 
after the spike. In the lower panel of figure [3] we present the measured SW flux curve with 
the dotted line and the modeled step function with the plain line. 

The resulting simulated lightcurves for O vn and O vm are shown in the upper panel of 
figure[3]in plain and dashed lines respectively. The plain vertical lines mark the limits of the 
total observing period and the dashed vertical line shows the limit between the ON and OFF 
cloud exposures. The dynamic model, taking into account SW measurements, is predicting 
a smooth but gradual rise of the oxygen line intensities in the ON and OFF exposures. The 
O vii and O vm line fluxes averaged for the ON-CLOUD exposure yield 3.56 LU and 0.5 
LU respectively, while averages on the OFF-CLOUD exposure yield 4.62 LU and 0.77 LU 
for O vii and O vm respectively. The equivalent averaged Suzaku data for O vn and O vm 
are shown in the figure with plain dots and hollow triangles respectively. 

The model-data comparison for the ON-cloud exposure shows that LB emission is neg- 
ligible, since the SWCX model accounts for all of the observed intensity. The X-ray flux 
increase (~30% and ~55% for O vn and O vm respectively) in the OFF-CLOU D exposure 
would change the data interpretation in terms of hot halo gas temperature, since lSmith et alj 
J2007h assumed a constant foreground emission. 
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Table 1 Comparison of shadowing observation foreground data to SWCX model for O VII and O VIII line 
intensities. " Observation date and observing instrument: C for Chandra, X for XMM-Newton and S for 
Suzaku. b MBM 12 Suzaku-off exposure is included in the analysis, as foreground e mission was not con- 
stant (s ee j|3.2.1l for details). c Revised values for SGF foreground line intensities, from Henley and Shelton 
f20073) . 







O VII (LU) 




O vm (LU) 




Target 


Observation " 


Data Data c 


SWCX 


Data 


Data 


SWCX 


MBM 12 
MBM 12 
MBM 12 
SGF 
SGF 


Aug. '00, C 
Feb. '06, S-on 
Feb. '06, S-off b 
May, '02, X 
Mar. '06, S 


1.79 ±0.55 
3.34 ±0.26 
5.68 ±0.59 

3.40 6.2+H 

0.13 


1.49 
3.56 
4.62 
3.16 
0.34 


2.34 ±0.36 
0.24 ±0.10 
1.01 ±0.26 

1.00 

N.C. 


l.Oil.l 


2.13 
0.50 
0.77 
1.02 
0.02 



3.3 Summary of Shadowing Observations 

In Table Q] we resu me the MBM 12 and Sou th Galactic Filament shadowing observations 
analysed in detail in lKoutroumpa et al. All four observations were simulated with the 

SWCX model with SW conditions as close as possible to measurements in solar instruments, 
to account for the SWCX contamination of the oxygen line emission in each field. 

Besides the basic differences due to solar maximum or solar minimum conditions dur- 
ing the observing period for each field, that we reported in section 13.11 three of the four 
observations (MBM 12 Chandra, Suzaku and SGF XMM observations respectively) were 
found to be highly contaminated by short-term SWCX emission variations. In particular, in 
our simulations we confirmed the contam i nation of the MBM 12 Chandra observation by a 
strong CME as suggested bv lSmith etai] J2005h . We suggested that the SGF XMM obser- 
vation was also highly cont aminated, most probably by a CME, which was confirmed by 
iHenlev and Sheltonl ( 12007 ah when they re-analysed and compared their Suzaku and XMM 
data. 

The observation data represent u nabsorbed (attributed to the LB) O vn and O vm line in- 
tensities as derived from the authors dSmith et allbOOibOOTlJHenlev et allboollHenlev and Sheltonl 



sis published by 


Henlev et alj 


Koutroumpa et al., 


2007 


)andt 



The shadows block more or less efficiently the Galactic Halo oxygen emission, and in 
their analysis the authors derive the unabsorbed foreground emission of the oxygen lines, 
that we use to compare to the SWCX heliosp heric emission. The only exception is in the 
Suzaku/MBM 12 observation dSmith et aflfeOOTt) where the O vn and O vm halo emission is 
added to the foreground emis sion for the OFF-C LOUD exposure. Contrary to the constant 
foreground assumption in the ISmith et al.1 d2007l) analysis, we demonstrated that there is a 
-30% and -55% increase in the OFF-CLOUD simulated SWCX O vn and O vn line inten- 
sities respectively, due to the brief SW enhancement at the end of the ON-CLOUD exposure 
(see 93.2. It . This increase was erroneously attributed to the Galactic Halo emission, which 
should be revised, and thus we include these values in the analysis as well. 
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O VII LINE ANALYSIS 

Old SGF analysis 

• data lit: 

a = (-0.26 ± 0.43) LU 
= 1.19 ± 0.14 

New SGF analysis 

o data tit: 

a = (0.09 ± 0.73) LU 
b = 1 .00 ± 0.21 



2 3 
SWCX Model Line Intensity (LU) 



O VIII LINE ANALYSIS 

Old SGF analysis 

A data fit: 

a = (-0.1 7 ±0.24) LU 
b = 1.19 ± 0.19 

New SGF analysis 

A data fit: 

a = (0.46 ± 0.44) LU 
b = 0.80 ± 0.38 



0.5 1.0 1.5 2.0 

SWCX Model Line Intensity (LU) 



Fig. 4 Linear fit s of oxygen line intensitie s. Foreground emission data from shadowing observ ations i n litter- 
ature (MBM 12: ISmith et ani200al2007l) and (SGF: lHenlev et aUl2007l;lHenlev and Sheltorll2007bl) versus 
SWCX model results for the O VII (left panel) and O VIII (right panel) line intensities. All linear fits are 
calculated with no error-bar weighting except for the O VII New SGF analysis fit (left panel, dashed line). 
See j|3.3l for more details. 



In figure[5]we linearly fit the data over the SWCX model results for the O vn (left panel) 
and O vm (right panel) lines for the initial (plain lines) and revised (dashed lines) data sets. 
We did not use data errorbars as standard deviation to weight the fit, because n ot all data 
errorbars were communicated (Henlev et all I2007I: IHenlev and Sheltonl l2007bh . The only 
exception is the O vn line fit, since Hen lev and Sheltonl d2007ah provided errorbars for both 
XMM and Suzaku O vn detections. The old and new fit coefficients (y(LU) = (a + bx)(LU)) 
for the O vn and O vm lines are also noted on the figure. 

The first linear fits (in plain lines in fig(4j we performed dKoutroumpa et all |2007|) had 
a slope of 1.19 (±0.14; 0.19 at lcr for O vn and O vm resp.) meaning that the SWCX model 
reproduced very well the local foreground emission measured in shadows. Moreover, the fit 
also gave an estimate of the residual foreground emission to be attributed to the LB. The 
foreground was found to be (-0.26 ± 0.43 at lcr) LU for O vn and (-0.17 ± 0.24 at lcr) 
LU for O vm which means that with a high probability the LB O vn and O vm emission is 
negligible compared to the heliospheric emission. 

The new fits (in dashed lines, figB>, where we substitute d the SGF observed foreground 
data with the revised values of IHenlev and Sheltonl d2007al) give slightly different results, 
although the correlation to the simulated values remains very satisfying (the correlation is 
even improved for the O vn line intensities) and the residual LB emission is still consistent 
to zero. 



4 SWCX in the 1/4 keV Band 

In this section we present a preliminary study of the SWCX background in the 1/4 keV (0.1- 
0.3 keV) band that we compare with the ROSAT and Wisconsin surveys data. In section |4~T| 
we describe briefly the SWCX simulation and in section l4~2l we compare the model results 
with the survey data sets. 
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Fig. 5 Calculated SWCX spectra in Line Units for the 1/4 keV band. Emitting ions are marked above the 
most intense lines. 



4.1 Simulations 



For the study on the 1/4 keV band, we updated our atomic database to include Fe, Si, S, Mg 
ions that emit intense lines in the 0.1-0.3 keV range. Individual radiative transition prob- 
abilities for these ions were calculated assuming that the ions are hydrogenic. Moreover, 
photon yields were calculated using a unique neutral species, which means that no dis- 
tinction between H and He was made. The hydrogenic ion assumption and unique neutral 
species choice are only approximations, but the results give a good estimate of the order of 
magnitude of the X-ray intensity for this energy range. An example of calculated spectra 
in the 0.1-0.3 keV range is presented in figure[5] with the emitting ion identifying the most 
intense lines. 

We have calculated spectra for two LOS, upwind (UW) and downwind (DW), for an 
observer assumed to be at 1 AU and in a crosswind position with respect to the LOS. This 
observing geometry was chosen to match the ROSAT satellite observing constraints. In the 
1/4 keV (0.1- 0.3 keV) range the model yields a total of 28.35 LU for the UW LOS and a 
total of 20.27 LU for the DW LOS. The intensity is probably underestimated because of the 
SW ion propagation in the heliosheath up to the heliopause, and in the heliotail up to ~3 000 
AU, where all ions are used up. The outer heliospheric region is neglected (for the moment) 
in our model, but preliminary estimates yield a maximum additional ~20% contribution in 
the DW direction. In general, because of parallax effects as well, we can assume that within 
a 50% error bar, these values give an good estimate of average values of the heliospheric 
SWCX emission in the 1/4 keV range. 
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SWCX model 


Data range 


Band 


UW 


DW 




B (cts s _1 ) 


11.99 


8.04 


[20 - 100] 


C (cts s _1 ) 


73.32 


56.24 


[50 - 250] 


Rl (RU ") 


141.14 


104.38 




R2(RU) 


192.07 


147.27 




R12 (RU) 


333.21 


251.65 


[250 - 820] 



Table 2 Summary of main band fluxes for the SWCX 
model in the upwind (UW) and downwind (DW) di- 
rections and measured data. For ROSAT data only the 
unabsorbed fraction I„ is given. " RU = 10~ 6 cts s 
arcmin~ 1 . 



4.2 Model-Data Comparison 

The sounding rocket survey of the University of Wisconsin and the ROSAT satellite survey 
have dominated the diffuse SXRB observations from 1972 to 1991 . The Wisconsin survey al- 



lowed a full-sky mapping of the IS SXRB with a spatial resolution of ~7° ( McCammon et all 
Il983h . The ROSAT survey improved the spatial resolution to ~12 ' dSnowden et alj,ll997h . 



However, the instruments did not have any spectral resolution so they provided only total 
fluxes for each band. 

In order to compare the modeled SWCX emission level to the measured SXRB inten- 
sities at 1/4 keV, we convolve the calculated spectra for the UW and DW LOS with each 
band response in this energy range of the instruments. The bands concerning our study are 
B (0.13-0.188 keV), C (0.16-0.284 keV) for the Wisconsin survey and Rl (0.08-0.284 
keV), R2 (0.14-0.284 keV ) for the ROSAT survey. The corresponding ban d responses are 
resen ted in figures: 1(a) of lMcCammon et all dl983h for Wisconsin and 1 of lSnowden et alj 



pit 
<2< 



2000) for ROSAT. In table [2] we summarize the total band fluxes derived from the SWCX 
spectra convolution with the band responses, and the measured data range for each band 
(evaluated from figures [6] and |7J. 

4.2.1 Wisconsin Sounding Rocket Survey 

The SXRB survey of the University of Wisconsin was performed with a series of 10 sound- 
ing rocket flights between 1972 and 1979. The only spectral information in the 1/4 keV range 
is extracted by the boron (B) and carbon (C) filt ers incorporated into the d etector windows 
that separate the two bands B and C presented in lMcCammon et all dl983|). figure 1(a). 



Detailed maps derived from the complete survey were presented in iMcCammon et al] 
dl983h . The maps present a clear negative correlation with the column density of galactic 
neutral hydrogen Nhi- T he X-ray intensi t y is low est toward the galactic plane and highest 
at high galactic latitudes. ISnowden et al.1 dl990bh investigated the band intensity ratio B/C 
(low to high energies) and found a variation between 0.25 to 0.46 in a dipole-like spatial 
correlation, aligned with a roughly galactic center-anticenter direction, that they attributed 
to an anisotropy of the Local Hot Bubble temperature ranging between 10 6 2 and 10 5 9 K 

respectively. 

From lSnowden et all dl990bl) we extracted a scatter plot (fig|6]> comparing the B and C 
band measurements in the Wisconsin survey. The data range is summarized in table |2] and 
gives ~[20-100] cts s for the B band and ~[50-250] cts s for the C band. We superposed, 
in this figure, the two points (wide circles) corresponding to our simulated spectra convolved 
with the B and C band responses (see table [2]l. Error bars are given only for information at 



12 



120 




50 100 150 200 250 300 



C BAND 



Fig. 6 Scatter plot of the total intensity calculated in the SWCX model (large circles) compared to total 
measured (crosses) values in B and C bands in the Wisconsin survey. The model reproduces the observed C 
band fluxes in the gala ctic plane, but not the B band fluxes. Figure extracted and adapted for our study from 
ISnowden etal] <1990bl) . 

50% of the SWCX band flux, in order to give an estimate of the heliospheric flux range in 
these bands. 

The SWCX model in the C band predicts fluxes well within the observed values in the 
galactic plane (lower limits), but simulated B band flux is about two times lower than the 
lower observed limits. The band ratio B/C for the SWCX flux is accordingly inconsistent 
with measured values, and fairly constant from UW to DW directions, with a mean value of 
B/C -0.15. 

4.2.2 The ROSAT Survey 

The ROSAT all-sky survey was performed in a scan mode with imaging telescopes in the 
soft X-ray band of 0.1 - 2 keV, during the first year of the mission. The observing geometry 
constraints imposed to scan the sky in a big circle, during one orbit, perpendicular to the 
spin axis that coincided with the Sun-Earth direction. Therefore a complete scan of the sky 
was performed in a six-month period. 

The ROSAT all-sky maps in the 1/4 keV range dSnowden etall 1 1 995l . [l997h confirmed 
the Wisconsin results and the negative correlation with the H I column density. The improved 
spatial resolution of ROSAT, with respect to previous surveys, allowed a better use of shad- 
ows as tools to separate the foreground from distant components of the SXRB. An extensive 
catalog of shado ws in the 1/4 keV ban d, with fitted foreground and distant emission values 
was presented in lSnowden et al. 

In figure [7] we present a scatter plot of foreground intensity I„ measured in band R12, 
which is the sum of intensities in bands Rl and R2. Units are in 10~ 6 cts s _1 arcmin , better 
known as ROSAT Units (RU).The figure was extracted from lSnowden et all feOOCh and the 
authors compare the data set analysed with two different methods. The exact comparison 
of the two methods is beyond the scope of this paper, since what we need is just the data 
intensity range (250-820 RU). In the figure y = x, thus the fitted slope is one. Simulated 
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Paper I Analysis Intensity 
(10~ 6 counts s _1 arcmin -1 ) 



Fig. 7 Same as fig.|6]but for the unabsorbed (I„) emission of ROSAT shadowing observations in band R12 
(Rl + R2). The model repro duces most of the R12 band flux in the galactic plane. Figure extracted and 
adapted for our study from Snowden et at (2000). 



spectra convolved with Rl and R2 band responses and summed for the total result in R12 
band yield 333.21 and 251.65 RU for the UW and DW direction respectively. These values, 
represented by the wide circles in fig. [7] with error bars at 50% (to give the SWCX flux 
range), are within the lower limits of observational data, that correspond to the low galactic 
latitude regions. 

In the same paper jSnowdenetailbOOfj) confirm the dipole-like spatial correlation of 
the R2/R1 (high to low energies) flux ratio, although the variation is less pronounced than 
the one derived from the Wisconsin analysis, 1.25 to 1.04 from UW to DW. Our model 
results yield a constant R2/R1 ratio at an average value of 1.38. 



5 Discussion 

We have presented simulation results on the SWCX soft X-ray emission in the heliosphere 
in the 0.1-1.0 keV energy range, that we compared to data from past and present observing 
instruments. We separated our analysis into two energy ranges: the 3/4 keV, dominated by 
the oxygen (O vn and O vm) lines and the 1/4 keV, where heavier ion lines as Mg, Si, S and 
Fe dominate the SWCX spectra. This analysis yields estimates of the SWCX heliospheric 
component within the diffuse Soft X-ray Background and confirms the large contamination 
of X-ray data by the heliospheric emission that needs to be taken into account in future 
determinations of the LB temperature and pressure. 

The simulation results on shadowing clouds observed with Chandra, XMM-Newton and 
Suzaku suggest that the local 3/4 keV emission detected in front of shadows is entirely 
explained by the heliospheric SWCX emission (see table Q] and figure [4]l and no emission 
from the LB is needed at these energies. Thus, combining our results in the 3/4 keV range 
and previous results published on the 3/4 and 1/4 keV ROSAT band ratios, the strict upper 
limit of the LB temperature is 10 6 K. Indeed, if the LB plasma is in collisional equilibrium 
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at 10 6 K, then it mainly emits in the 1/4 keV range, and only very little emission is produced 
in the oxygen lines at 3/4 keV. 

But, then, how does SWCX affect the data interpretation in the 1/4 keV range? We have 
compared SWCX simulation results with data from the ROSAT and Wisconsin surveys. 
Both sets of measured data were originally interpreted in terms of hot gas emission in the 
LB at a temperature of ~ 10 6 K. Variations in the observed X-ray intensity were attributed 
to variations in the extent of the emission volume and therefore the emission measure of the 
plasma. 

SWCX emission flux calculated for the R12 and C bands (which practically cover the 
same energy range) yielded average values that explain most of the observed emission in the 
galactic plane, previo usly attributed to the restricted regions of the LB. Indeed, according to 
ISnowdenetaT] ( ll998h . in the galactic center and anticenter directions (roughly UW and DW 
respectivel y), the LHB is found t o have an extent of Rlhb ~ 50 pc (see upper left panel of 
figure 10 in lSnowden et all 1 19981) . which corresponds to an unabsorbed 1/4 keV emission of 
I = Rlhb / 0.155 » 322 RU measured in the R12 band, while our model predicts ~[250-333] 
RU. 

On the other hand, the SXRB intensity is typically brighter at higher Galactic latitudes, 
and SWCX cannot account for all the emission. This may be explained by the fact that at 
high latitudes the LOS's point through the chimneys allowing for non-absorbed halo emis- 
sion to reach the observer. Moreover in the B band, which is the least overlapping in energy 
range with respect to the other 1/4 keV bands, the SWCX simulations predict about half of 
the observed lower limits. Besides, the SWCX model fails to reproduce the observed R2/R1 
ratio by 10-33% and the B/C ratio by 40-67% . 

The answer, then, to the local unabsorbed SXRB puzzle might lie in the proper mix of 
SWCX emission and warm (rather than hot) LB that will fill-in the gap in the B band and 
reproduce the observed band ratios (Koutroumpa et al., in preparation). 
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